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Dissipative effects on the nonadiabatic transition for the two and three level systems are studied. 
When the system is affected by a strong dissipation through the diabatic states, the exact transition 
probability is enumerated making use of the effective master equation. In the two-level system, we 
consider the case where the external field is swept from not only a negative large value but also from 
the resonant field, and the exact transition probabilities in these cases are derived. The transition 
probabilities are derived for the three-level system where the three diabatic states form only one 
avoided level crossing point. These probabilities are compared with the one in the pure quantum 
case obtained by Carroll and Hioe. 
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I. INTRODUCTION 



Nonadiabatic transition at an avoided level crossing point plays a crucial role in quantum dynamical changes of 
states, and yields the variety of phenomena in physics and chemistry. The well-known Landau and Zener (LZ) transi- 
tion probability clarifies the roles of the energy gap and the sweeping velocity of the external field in the nonadiabatic 
transition of two-level system Although the LZ transition probability is given in the two- level system, it is 

approximately applicable to multi-level system where the avoided level crossings are effectively well described by only 
localized two levels. Hence it is adopted in the analyses of many experiments which treat time dependent phenomena, 
such as collision of particles (||| i optics , and magnetic phenomena [^jl3| ■ For general multi-level systems where 
many levels can simultaneously affect on each other, different formulas of transition probabilities are derived for sev- 
eral models, i.e., the one where only one level interacts with a band of levels |TJ-|l^], the generalized model for this 
model (f7]jli| and the bow-tie model where many levels form only one avoided level crossing [|l9|j2C| 1 . We should also 
note Brundbler and Elser's hypothesis which states that the survival probability of the diabatic state with maximum 
or minimum slope is described by the exponential form determined by only the velocity and the off-diagonal elements 
in the Hamiltonian [ pf[ . 

On the other hand, we must also consider the effect of dissipation, since real experiments are always exposed 
to thermal environment. The thermal environment causes decoherence and the inevitable deviation of transition 
probability from the one of pure quantum case occurs. This modification becomes significant in real experiments such 
as the adiabatic rapid passage with phonon couplings p2| , the nonadiabatic transitions in localized centers in solids 
PH , and nonadiabatic magnetization process in molecular magnets such as Mni 2 and Feg |24|-|27|] . Kayanuma studied 
such thermal noise effect for the two-level Landau-Zener model, and derived a formula for the effective transition 
probability in the limiting case of strong damping dissipation by perturbative approach p8|] . The effective transition 
probability becomes 1/2 in the adiabatic limit due to the dissipation effect, whereas it converges to the asymptotic 
expression of the LZ probability in the fast sweeping limit. Ao and Rammer carried out first principle calculation 
to investigate temperature dependence of the transition probability of the two-level system with phonon reservoir 
which corresponds to the Kayanuma's situation in the high temperature limit in case of the Ohmic spectral density. 
Especially they found some compensation effect that the transition probability for zero temperature becomes the same 
value as the Landau-Zener probability p9| . 

In this paper, we study such thermal noise effect in not only the two-level system but also the three-level system. 
Thereby we try to investigate the effect of multi-level with thermal noise. We exploit the method to analyze strong 
dissipation effect using an effective master equation instead of the perturbation approach adopted in previous studies 
p3| , p8|j29| . We show that the effective master equation approach is very convenient for deriving the transition proba- 
bility in the strong damping limit. Using this approach we first reproduce Kayanuma's formula in the two-level system 
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when the external field is reversed from large negative value to large positive value. We next consider the situation 
where the field is swept from resonant field (zero field) to large positive field, and derive exact transition probability. 
As a result the exact relations between these cases are found. The three-level model we consider is the same model 
as Carroll and Hioe considered In the model, three diabatic states form only one avoided level crossing point. 
Therefore the transition mechanism is quite different from the LZ mechanism which describes transitions between 
local two levels. Therefore we see the effect of multi-level not only in pure quantum case but in dissipative case. We 
adopt the effective master equation approach and derive the transition probabilities in the strong damping limit. The 
probabilities are compared with the one of pure quantum case obtained by Carroll and Hioe. The expression of the 
probabilities are always the same regardless of level structures, although in pure quantum case the expressions of the 
transition probabilities show some variations. 

This paper is organized in the following way. In section II, we derive the master equation when strong noise couples 
with the system and compare it with the master equation for the system with phonon reservoir. Section III is devoted 
to the problem for two-level system, and we derive the transition probabilities for three-level system in section IV. 
Summary and brief discussion are given in section V. 



II. MASTER EQUATION 



We derive the master equation for the system with dissipation. Throughout this paper, we study the transition 
probability by solving the master equation. The master equation we shall consider is derived for various types of 
dissipative environments such as stochastic noise field and phonon reservoir. We here choose a stochastic noise field 
as a source of dissipation and rigorously derive the master equation. As shown in Appendix A, this master equation 
can be obtained in case of phonon reservoir with the Ohmic type spectral density at very high temperatures. The 
correlation of the stochastic noise is assumed to be very short. This situation is described as 
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(1) 
(2) 



where is a noise which affects on the system through the £th operator Xg. The matrix Xi is diagonal in the 
diabatic bases of the Hamiltonian TL(t) so that the computability [Xi, Xi/\ = is satisfied for arbitrary I and £'. The 
noise is supposed to be the white Gaussian process. We start with the Von-Neumann equation for the density matrix 
in the interaction picture (h = 1 here and hereafter), 



dt 



$>(*)A(t)p (I) 



(3) 

(4) 
(5) 

(6) 

Here exp<_ and exp^ express the time ordered product of exp onentials. In the case of white Gaussian process (^), 
there are several approaches to derive master equation |30| , |3l|| . Here we use the Novikov's relation, which holds for 
arbitrary function <?([£],£) | pH , 



p®(t) = cxp^_ ^— i J duH(u)^ p(t)cxp^ J duH(u) 
Xf (t) = exp^_ ^— i J duH(u)^ Xgcxp^ J duTL{u) 



(7) 



where the symbol [£] means that <?([£], t) is a function of the process of noise £, and (...) means the average over the 
noise £(t). By use of this mathematical formula, the average over noise for Eq.(0) is reduced to, 



d(p {1) (t)) 
dt 



*(«) 
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£ Jt ° \ t< Jt ' J 

= J2^i(t)iP m (t))- (8) 

I 

Here we used Novikov's relation and the properties of the white noise. Thus we arrived at the master equation in the 
Shrodinger picture of the density matrix 

^ = -i [H(t),p(t)} [X t , [X e , p(t)}] , (9) 

in which we denoted p(t) for (p(t)), omitting the symbol of average over noise ( ). In the following sections, we 
focus on the noise with short correlation (^|) and large amplitude of 7, namely the strong damping (SD) limit. This 
is equivalent to the phonon reservoir with high temperature. We investigate the properties of the nonadiabatic 
transitions based on the master equation (||). 

III. DISSIPATIVE EFFECTS IN TWO-LEVEL SYSTEM 
A. Transition Probability 

In this section, we consider the noise effect in the two-level system. For two-level system, the transition probabilities 
in the strong damping case have been studied by several authors using perturbation approach |^,^]. Here we 
alternatively adopt the different approach, deriving the effective master equation which simplifies derivations of 
transition probabilities for the case. Here two situations are studied. We first consider the familiar situation that 
the external field is reversed from —00 to 00. In this case, we demonstrate that Kayanuma's transition probability 
Pgo" is reproduced easily using the effective master equation. We second consider somewhat unfamiliar but realizable 
situation that the external field is swept from to 00 when initially a diabatic state is occupied. In magnetic system, 
this situation corresponds to the case where a spin is saturated to the down (or up) state under very strong field and 
next the field is switched off and the field is swept linearly in time from zero field. 

The model we shall consider in this section is described as, 

H t at(t)=H(t)+Z(t)a z , (10) 
H{t)=-vt\o* + T 1 -a*, (11) 

where a a is the a(— x, y, z) component of the Pauli matrix. The diabatic states correspond to the down state |1) and 
up state 1 2) which satisfy <r z \l) — — 11) and a z \2) = |2), respectively. T is the transverse field which is responsible for 
the tunneling between the diabatic states. Here we take only a z as the operator on which the noise acts in (j^). Thus 
the master equation (||) is concretely written as, 

§- t P(t) = ~i\ [vta z + Ta x ,p(t)} - I (p(t) - a zP (t)a z ) . (12) 



We define the following variables: 



ci = pu - P22, (13) 

C2=Pl2, (14) 

c-3 = P2i- (15) 



The time evolutions of these variables are determined by the differential equations: 



ci = -ir(c 3 - c 2 ), (16) 
iT 

C2 = (-ivt - 7)c 2 + —ci, (17) 
iT 

c 3 = {ivt - 7)c 3 - — ci. (18) 



3 



Here we consider the SD limit, 7 — > 00. We first consider the variable c 2 (t), which is formally solved from Eq.([L7| 
We can approximate it by partial integral; 



02 (*) = c 2 (t ) + l -L e - mt2 / 2 ~-' t I due lvuA l 2+ ^ u Cl {u) 



In 



c 2 (t ) + %e-^ 2 /^ 




vu /2-\--yu 



c\{u) 
ivu + 7 

d_ ( ci(t) \ 
dt \ ivt+~{ J 



I 



, ivu 



72+711 " 



<i / c\{u) 



^ 2v J t du \ivu + y 



ivt + 7 ivt + 7 



. . iT ci(t) , „. 

Here we used the fact that the term e -7 ^"* -' is negligible due to large 7, and we neglected the higher order terms of 
(ivt + 7) _1 [^2). When the diabatic state |1) is initially occupied, namely, 

c l (*o) = l, c 2 (t ) =c 3 (to) = 0, (20) 

C2(t) and C3(t) are approximated as, 

ir ci(t) 

C2W = TiStT? (21) 

c 3 ft) = — Cl(t) . (22) 
2 — ivt + 7 

These relations lead us to the simplified equation for the diabatic states. By substituting the relations (|21|) and (|2 
into the equation (|l6|), we arrive at the effective master equation : 



2iw [i + 17/11 t — vy/v 

Now let us consider the first problem, i.e. to = —00. In this case, we can readily integrate the master equation (|2 
to get, 

ci(oo) =exp • ( 24 ) 

We now consider the tunneling probability from the state |1) at t = —00 to the state |2) at t = 00. This corresponds 
to the value of /c»22(oo)- By using the conservation of probability Trp = 1, this transition probability PJ^ (= ^22(00)) 
is obtained, 

^-5(1— p (-£))■ <*> 

This is nothing but Kayanuma's transition probability psfl 

For the second case where the field is swept from resonant field, i.e., to = 0, ci(oo) is readily calculated as 

ci(oo) =exp (-^— Y (26) 
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which gives the transition probability , 



The validity of these probabilities ( |25| ) and (27) is confirmed by numerically integrating the master equation (|12|). In 
figure 1, we present numerical data and the analytical results given in ( |25| ) and (^) as a function of sweeping velocity. 
Here V is 0.01, and 7 is 10.0. We show not only the dissipative case but also pure quantum case (P{^~ and )• We 
see that formulas (p5|) and d27]) agree with the numerical results almost perfectly. 



4 



1 


Plz + 

^SD x 


1 

-\ 


Plz 

PSD ° 


h \ 




® k \ 
~ \ \ \ 








" '^\%^ 









1 23456789 10 11 



FIG. 1. Transition probabilities as a function of sweeping velocity v. Parameters are taken as T = 0.01, and 7 = 10.0. Points 
are numerical data, and lines are the analytical values. Here the subscript 'LZ' means that it is pure quantum case. P°z ^ s the 
probability in the case where the external field is swept from zero field in pure quantum case. 

The variable c\ (t) is directly connected with the magnetization M(t) — Trcr z p(t). We obtain the magnetization 
process solving (p3|) as, 





■(*) = 


— exp 


T 2 

V 


M 0H 


■(*) = 


— exp 


T 2 

V 



L ( 7 \ IT 

— arctan — 

v \vtJ 2 



(28) 
(29) 



where the function y = arctanx is defined in the region of x € [— 00,00] and y € [— 7r, 0]. This shows that the the 
magnetization process depends on the noise strength 7, whereas final magnetization does not. This is also numerically 
confirmed. 



B. Relation between Psd and Plz 



The exact relation of the transition probabilities for pure quantum case and in strong damping case is discussed. 
When the field is reversed from large negative value without dissipation, the transition probability is denoted by ■ 
Here the subscript 'LZ' means that it is the pure quantum case. In the case of the field swept from the resonant point 
without dissipation, the transition probability is written as P^z • As briefly shown in Appendix B, the transition 
probability P°z 1S calculated as, 

« -I ('-« » (-£))■ < 3 "' 

Thus transition probability in each case is written as follows, 



P-+ - I 

^SD — 2 



1- exp P£-+ = l-exp(-^) 



P°+ = I (l - exp (-£)) P°+ = I (l - exp 
These asymptotically exact probabilities satisfy the following relations (see also the Fig.l) 



( 31 ) 



^sd + < ^lz + , (32) 
> Pit (33) 



The inequality ( |32] ) means that the dissipation reduces tunneling that the state remains in the ground state. When 
v -C r 2 in the pure quantum case, almost adiabatic tunneling from |1) to |2) occurs, whereas in the presence of 
dissipation, thermal excitation from ground state to the excited state represses such adiabatic evolution. Thus the 
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inequality ([52]) is easily understood. On the other hand, the inequality (|3^) ibdicates the opposite property. This 
effect is intuitively explained as follows. The initial state ip(0) (= |1)) is the superposition between the ground state 
|G(0)) and the excited state |.E7(0)) , 

m = \(\G(0)) + \E(0))) (34) 

In case of almost adiabatic evolution in the pure quantum case, v <C 1, the state of the system almost follows such 
superposition at t, 

ip(t) ~ i (eW*\G(t)) +e* (t) |S(t)>) . (35) 

where e* (t) and e* (t) are the dynamical phases. Since \G(t)) -> |2), \E(t)) -> |1) in the limit of t -> oo, the 
maximum value of transition probability is \ in the pure quantum evolution. In the presence of the dissipation, the 
noise also induces such uniform distribution, because the dissipation we now consider can be regarded as thermal 
effects with very high temperature. As a result the tunneling probability is larger in presence of the dissipation. 
We may say that the inequality (|3^) is a consequence of the special initial state, because if the initial state is not a 
diabatic state, e.g. ip(0) = |G(0)), the inequality ([33]) is not realized. We expect these characteristic relations ( |3^ ) 
and (|3j|) will be verified in real experiments. 



IV. THREE-LEVEL SYSTEM 



We here apply our approach to the multi-level systems, and investigate the dissipative effect on many levels. For the 
sake of definiteness, we focus our attention on the properties of the three-level system, where all diabatic levels form 
only one avoided level crossing point. Because the avoided crossing structure is not formed by localized two states, 
the transition mechanism is very different from the LZ mechanism. Therefore this three-level system will provide 
much information about the effects of strong correlation of many levels. This model was first studied by Carroll and 
Hioe JL9[ , and the exact transition probabilites have been obtained. The formulas for the probabilities show some 
varieties of expressions according to the relation between the slopes of diabatic states. The generalized arbitrary 
N- level system for this model is called the bow-tie model pPf ], and the transition probabilities and characteristic 
mechanisms of transitions are discussed |B3| . Thus this model is quite convenient for comparison of the transition 
probabilities in the pure quantum case and the dissipative case. We should also note that there exist some proposals 
for physical realization of this model by using optical systems • 



A. Analysis of Transition Probabilities 



We here enumerate exact transition probabilities in the three-level system in the presence of dissipation. The 
Hamiltonian we consider is written for the diabatic bases: 

H(t)=Ho(t)+Y,Zk(t)X k (36) 

(a™ \ 
X k = 4 fc) (37) 
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Here the kih white noise acts on the fcth operator X k . In this case, the matrix element of the master equation ||) is 
written in the following form, 

PI T} - = -i {H Q (t)e k p k m (t) - pe,k(t)Ho(t)k,m) ~ %mPe.,m[t), (38) 
at 

where 7i jm is written using the matrix elements of the operator X k and the amplitude of the noise j k as, 
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-«{ r i(Pio 


- poi) ^ 


- 1^2 (P20 


- P02)} , 




pll(t) = 


-*ri(poi - 


Pw), 










-«r 2 (po2 - 


Pio), 








Poi(t) = 


-i{Ti(pu 


- poo) ^ 


" T2P21 


- fiipoi} 


- 701P01, 




-i{r 2 (p 2 2 


- Poo) ^ 


" Tipi2 


~ V 2 tp02} 


— 702P02, 


M*) = 


-i{TiP02 - 


- r 2 pi - 


f («1 " 


V 2 )tp 12 } - 


- 712P12- 



Thus the differential equations for all matrix elements are given by, 

(40) 

(41) 
(42) 
(43) 
(44) 
(45) 
(46) 

Here we confine ourselves to the SD limit 7^ — + 00, and the initial condition poi(^o) = Po 2 (to) = Pi2(*o) = 0. In the 
same manner as in the two- level case, pi 2 (t) is approximated, expanding by partial integrals, 

p 12 (t) = exp j u t j J du exp ^ h 7 12 u j [iT 2 p w (u) - iT 1 p 02 {u)\ 

_ iT 2 pw(t) - iT 1 p 02 (t) ^ 
i(vi - v 2 )t + 712 

Po 2 (t) and poi(i) are written as, 

P02W = exp 702^ rfu exp 

-^2 (P22^) - Pooft)) - gVi2_ft) 

-w 2 t + 702 

Poi(i) = exp 7oi*J du exp 

^ -i^i (Pnft) - Pooft)) - iT 2 p 21 {t) 
-ivit + 701 



-i—z 1- 702^ I {-iT 2 (p 22 (u) - poo(u)) - iFipi 2 (u)} 



-i— 1- 7oim ) {-iTx (pii(w) - Poo(")) - «r 2 p2i(M)} 



(48) 



(49) 



The terms p\ 2 {t) and p 2 \(t) in (|4q ) and (49) are negligible because of Eq.([47|). Therefore we can approximate ( |4S| ) 
and d49|) as, 
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v 2 J 
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Vl ) 




^ jToi 

Vl 



ml 1 • = ( — )( —- I (paa(*) - Poo(*)) , (50) 

5ii(t)-poo(*)). (51) 



Thus we arrive at the effective master equations defining ci(t) = pn(i) — poo(t) and c 2 (t) = p 22 (t) — poo(^) : 

-2*r? (1 1 \ »r| / 1 1 \ 

ci = 1 : 1 - / : — ^-j— c i " 1 , - / - : — ^-7— c 2 ( 52 ) 



ui V t + *7oi /Vl t — ijoi/vi J v 2 \ t + 1702 /v 2 t — 1702 / v 2 

-ir? / 1 1 \ 2ir| / 1 1 

^2 = 1 T I - / T ^"7" Ci M - ; ._ , ^-p- C 2 . (53) 



«i \t + i7oi/«i t — i'yoi/vij v 2 \t + i^ 02 /v 2 t-i^ 02 /v 2 

All matrix elements Pk,i(t) are obtained if poo(i), Pn(i), and p 22 (t) are calculated as solutions of these effective master 
equation. 

Now let us solve these differential equations. We start with the parameters which satisfy the special relation, 

701 702 

— = — =a. (54) 

M W 

Here a is always positive since 701 > and 702 > 0. In this special case, the equations are simplified in the form, 

^W ci $y (55) 



dt V c 2W y V * + j-W I c 2 (t) y ' 
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where the matrix M is given by 



M = 



2a b 
a 2b 



T 2 T 2 . 
a = - — - and b = -. — r 

M \v 2 \ 



(56) 



The right hand side of eq.(^5() has the explicit time dependence only in the prefacer. Therefore by diagonalizing the 
matrix M, we can obtain the scattering matrix which connects c(— oo) with c(oo). The matrix M has these eigenvalues 



A± = a + b ± \Ja 2 + b 2 - ab. (57) 
Using a, b, and these eigenvalues A±, the final state and initial state are connected using the scattering matrix S as, 



ci(oo) 
c 2 (oo) 



5*1,2 
"5*2,1 
5*2,2 



o | ci(-oo) 
1 c 2 (-oo) J ' 

1 



2Va 2 


+ b 2 
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— ab 


2Va 2 


+■ b 2 
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— ab 


2Va 2 


+ b 2 
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— ab 


2Va 2 


+ b 2 


— ab 



(a - b)(e- 27rX + - e- 2lrX -) + V a 2 + b 2 - a6( £ - 27rA + + e" 2 ^- 

2ttA_ 



( e - 2TA + -er 2lxX -) 



(a - b){-e- 27TX + + e- 27!X -) + V 'a 2 + b 2 - ab( e - 27rX + + e~ 27rA -) 



(58) 
(59) 

(60) 
(61) 
(62) 



Here a does not appear because it only gives the singular point in the Cauchy's integral to yield (p9[)-(p2[) , which is 
the same situation as in the two-level system. Therefore the scattering matrix does not depend on the concrete values 
of dissipation strength, 701 and 702 as far as those are large and the relation ( pi] ) is satisfied. We now obtain the 
probabilities for various initial states. 
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P SD + (0^0) 
^sd + (0-1) = j P S d + (1-0) 
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(63) 
(64) 
(65) 
(66) 
(67) 
(68) 
(69) 



In the adiabatic limit v\ — > +0 and v-i — > +0, all probabilities converge to 1/3 due to the strong effect of dissipation. 
As shown numerically in the next section, these formulas are always valid even when the relation (|54|) is not satisfied. 
That is, the probabilities are little affected by the variation of the strength of dissipation 701,7012, and 712- 



B. Numerical Investigation 



We numerically integrate the equation (p8|), and compare with the asymptotically exact transition probabilities 
obtained above for various parameter values. We write the slopes of the diabatic states v\ and V2 using parameter v, 



vi = a\v 1 
v 2 = a 2 v. 



(70) 
(71) 
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Dimcnsionless parameters oi and 02 gives the ratio of v\ to V2 and determine the level structure. We consider three 
types of level structures, namley (a) ax ■ a.i < 0, (b) ot\ > 012 > 0, and (c) «2 > ot\ > 0. In Fig. 2, the transition 
probabilities Psd~(0 — > j) (j = 0, 1, 2) are shown for these cases. The probabilities are plotted as a function of the 
parameter v for the parameters Ti = 0.1, T2 = 0.2, and various sets of (701,702,712)- Here the sets of (0:1,0:2) are 
taken as (1,-0.5) for Fig. 2(a), (1,0.5) for Fig. 2(b), and (0.5,1) for Fig. 2(c), respectively. The lines are theoretical 
values for the SD limit (0 — > j) given by (|6^ )-(|69|) and the probabilities for pure quantum case -Plz + (0 — * j) which 
are already obtained by Carroll and Hioe fll9[ ] . The analytical solutions of probabilities in the pure quantum case are 
listed in Table I. As can be seen in these figures, we can see good agreements between the numerical data and theories 
(|63|)-(|69|). We find little dependence on the variety of (701,702,712), that is, the formula (|63|)-(|69|) are valid even if 
the relation (|54|) is not satisfied as far as the dissipation is very strong. This was also confirmed for -Pgr/ (1 — > j), and 
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FIG. 2. Comparison the numerical calculation with theories for various sets of {7^}. Data are plotted as a function of v. 



Points are numerical data. 
Fig.2(b): (1,0.5), Fig.2(c): 



and the lines of P sr y and -P LZ ( see Table I) are theoretical values. Fig. 1(a): (0:1,02) = (1, —0.5), 
(0.5,1). 



As seen in Table I, the analytical expressions in pure quantum case f£z (0 — * j) show some variations according to 
relations of V\ and V2- On the other hand, the probabilities in the dissipative case (|63|)-(69) do not depend on such 
level structures. For instance, in cases of Fig. 2(a) and Fig. 2 (b) where the sweeping velocities are (v%, V2) = (v, — 0.5i>) 
and (1)1,1)2) = (v,0.5v), respectively, each probability P^~(i — * j) for the both cases completely agrees with each 



other because formulae ( 63 )-(p9|) have the same values for different (1)1,1)2) with the same absolute values. However 
in the pure quantum case, the probabilities -P£z + (* —* j) are different between these cases as found in Table I. This is 
a remarkable contrast between the dissipative case and the pure quantum case. 

In the slow sweeping v <C 1, the deviation of P^ from P^ are large. In the fast sweeping region v ^> 1, P^ 
asymptotically converges to the behavior of P{^~ ■ This means that the system is little affected from dissipation for 
fast sweeping because the time that the system stays at around avoided level crossing point is very short. This is the 







same behavior as found in the two-level system in pA 



TABLE I. The transition probabilities for quantum case. Here P = exp ( — nTi/vij and Q = exp f— tiT 2 ./^)- 





P c 7 + (0 -> 0) 


PJ+(0^1) 


P^+(0->2) 




PJ + (1^2) 


Pc + (2-2) 


Vl ■ V2 < 


(1-P-Q) a 


(1-P) (P + Q) 


(l-Q) (P + 0) 




(l-p) (l-Q) 


o 2 


\Vl\ > \V2\,Vl ■ V2 > 


P 2 Q 2 


(1-P)(1 + PQ) 


P(1-Q)(1 + P0) 




P(l-P)(l-Q) 


(1-P + PQ) 2 


U2| > • «2 > 


P 2 Q 2 


Q(1-P)(1 + PQ) 


(1-Q)(1 + PQ) 


(l-Q + PQY 


Q(l-P)(l-Q) 


o 2 



V. SUMMARY 



In two and three-level system, we derived the effective master equations which well describe time evolution of 
the system in the SD limit. Thereby we obtained analytical transition probabilities. The effective master equation 
approach is quite useful because the differential equation of system's variable becomes very simple. This approach 
will be applicable in the other systems whose exact transition probabilities can be analytically enumerated in the pure 
quantum case. 

In the two-level system, we consider the two cases where the external field is swept from large negative field and 
from zero field. Both situations are easily realized in real experiments. We hope that the exact relation (|3|) and (|33] ) 
is confirmed in real experiments using classical optical system and Cooper pair f35fl , and so on. 

The transition time in the two- level systems has been discussed in the literatures ]28|j3rj^3|] when the external field 
is reversed from large negative field. According to Vitanov's definition [Q, the transition times t tT is written as, 

tr = P22M _ Ci(oo) 

" P ' 22 (0) ~ c[(0) ' 172 j 

under the initial condition of pn(— 00) = 1. Vitanov derived the exact expression of transition time i^'z m the pure 
quantum case as follows, 



_ Vl - e-^/v 
LZ (I»cos( X ) ' 



(73) 




X= I + -gr(---)-argr^l--J, (74) 
where T(x) is the Gamma function. Eq.([73|) converges to 

r 2 „ 1 

(75) 

This asymptotic behavior agrees with the one in the earlier studies |28| , |36| . On the other hand, the transition time of 
dissipative case ig r D is readily written for the definition (|72|); 

T 2 / 71T 2 \ 

4 r B = T e X p(-— ). (76) 

This means that the transition time decreases as 0(1/7) with the increase of 7. In the three-level system, the 
analytical expressions (|6^)-([39|) are the solutions under the condition (|54|). However these solutions are valid beyond 
the condition (^34]) as shown numerical calculation. According to Carroll and Hioe's analytical solution (Table I), 
the pure quantum transition probabilities show some variations depending on the level structures. However, the 
probabilities in the dissipative case do not show such dependences. It is interesting to confirm experimentally this 
thermal effect because the three- level system we consider can be realized experimentally fll9| , |34"| | . 

Eq.(||) is derived under the condition where the noise affects on the system only through the diabatic states. In 
case of the off-diagonal coupling, i.e., Xi = a x , we also derived the similar master equation. In this case, we can easily 
show that in the strong damping limit, all the transition probabilities become uniform regardless of initial condition 
p5| , ^8[ . Thus the transition probability is affected by the coupling form with the thermal environment. Therefore it 
is also interesting to consider the transition probabilities for various coupling forms with finite 7. 
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APPENDIX A 

The same type of master equation as Eq. ( p5| ) is derived in the case where the system of interest couples with the 
phonon bath through the diabatic states as follows |39|, 

«« = WW + A £ X t £ la (ftCf)t + 6 W) + J2 ub^bW , (77) 



where H(t) denotes the system Hamiltonian of interest, and Xg is the £th operator which interacts with the phonon 
system. We assume the computability between the coupling operators Xg, i.e, [Xp,X^] = 0. The operator b$ 

(£)..... 

and bcj is the creation and annihilation phonon operator which interacts with the system through the £th coupling 
operator Xg. 

We use the projection operator technique to trace out the reservoir's degree of freedom ]4C|]and assume that the 
correlation between reservoir's variables is short-lived. Then we obtain the master equation for the system in the 
second order of coupling strength A pl[^2[ ; 



d 



p(t) = ±[H, P (t)]-\ 2 J2 T tP( t )> (78) 



dt ih 

t 

where Tgp(t) is given by 

Fep(t) = -2 dt' due iut '$ e (u){X e X e {-t , )p(t) 
n Jo J-oo 

- ef**"X iP (t)X t (-l!) + e^p(t)X e (-t')X e - X e (-t')p(t)X e ) , (79) 
Here Xg(~t') means the Heisenberg operator at time — t', 

X e (-t') = exp^ {^\J duH{u)\ Xexp^f- J duH{u)^j (80) 



In case of phonon reservoir described in (|77|), is given by, 

*/H = R — ePhL0 _ 1 — , (81) 

where (3 is an inverse temperature 1/T. Ii(lo) is called the spectral density. 
We restrict ourselves to the case of the Ohmic spectrum, 

h(u)=hw, (82) 

and high temperature, 

T > 1. (83) 

In this case, by using the fact, 

-» HI e T, (84) 



the master equation (78) is reduced to the simple form, 

- ^ [H,p(t)] - \ 2 TJ2^ [Xt, [X e ,p(t)}} . (85) 

Although the equation (|7^) derived by the projection operator approach is an approximation because higher order 
terms of A is neglected and the fast relaxation of reservoir is assumed to make the equation Markovian, it can well 
describe features of evolution of the system especially in case of high temperature. 
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APPENDIX B 



In this Appendix, the transition probabilities ( |3l"| ) are derived. We start with the Hamiltonian, 



H(t) = - {Ta x - vta z ) . 



We consider the Schrodinger equation. 



ih^(t) = - (Ta x - vta z ) *(<). 
dt w 2 v ' x ' 

This equation ( |87| ) is concretely written defining the component like ^(t) — (x±(t),x 2 (t))^ , 

vt r 

i%i(t) = —xi(t) + -x 2 (t) 

r vt 

ix%{t) = -xx(t) - — x 2 (t). 



Now we transform the variables to the following ones, 



xi(t) 



then we obtain the equation, 



2/2 (t) := x 2 (t), 



2 lT - yi + {l -- )yi --y 2 =0 

n .d v r 

2t-y 2 + -2/ 2 - -2/x = 0. 



Here we used the relation d/dt = 2td/dr. By using the Laplace transformation |44|] defined as, 



we obtain the integral representation for x\ (t) and x 2 (t) after straight- forward calculation, 



2 ~* 8v 



9 1 <^r HiSI 



Here the integral counter must satisfies the following condition, 

"(2ic-|)i,(0^ T 

For the variables £ = |£|e l( ^, t = |£|e* e , we choose the counter with the condition, 

6 + 26 = 3tt, 



c e 



(86) 
(87) 

(88) 
(89) 

(90) 

(91) 
(92) 

(93) 
(94) 

(95) 

(96) 
(97) 

(98) 
(99) 



for large |£| noting the relation £i = |£||t|e l ^ +2e ). When the initial time is t = oo • e l7r and the final time is t = oo- e*°, 
the phase of <f> varies from it to 3tt from the relation (|99|). We now consider the initial condition as, 



xi(— oo) = 1. 



(100) 



This condition is realized in the contour which encircles the singular point £ = — ij and choosing the constant A as, 
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r(0+) 

A = e^/ I du(-u)-^ 2+ir ^ Sv e- u . (101) 



Thus the wave function at t — oo is calculated by using analytical continuation following the condition (|99|), i.e. 
7r — * <j) — * 37T Q. Thus revival probability is calculated as, 

^(|t|e i0 )=^(|t| e iw )e X p^-^ (|t| - oo) (102) 

This means nothing but the relation of Landau-Zcncr transition. 

Next we consider the case where the external field is swept from zero value [[l9],[[^, that is, 6 = and +0 — > |t| — > oo. 
In order to derive the survival probability in this case, we first note the time symmetry that the probability is the 
same value as the one obtained when the initial time is taken as t — — oo (the external value = — oo) and the final 
time is t — —0 (the external value = —0). We consider the latter case (— oo — » t — ► —0) because we can use the 
same contour as the previous case (|9^) which satisfies the initial condition ( |lOC| ). We can readily write the integral 
representation of the x\{— 0) and obtain the transition probability P^, + : 



(-0) = A\t\ (-^-WrVs* (_ x _ ^ lT2/Sv e — W4|t|3 



.4 



1 + e 7rF2 /' kv 
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